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Abstract. Random walks on complex networks, especially scale-free networks, have attracted considerable 
interest in the past few years. A lot of previous work showed that the average receiving time (ART), i.e., the 
average of mean first-passage time (MFPT) for random walks to a given hub node (node with maximum 
degree) averaged over all starting points in scale-free small-world networks exhibits a sublinear or linear 
dependence on network order N (number of nodes), which indicates that hub nodes are very efficient in 
receiving information if one looks upon the random walker as an information messenger. Thus far, the 
efficiency of a hub node sending information on scale-free small-world networks has not been addressed 
yet. In this paper, we study random walks on the class of Koch networks with scale-free behavior and 
small- world effect. We derive some basic properties for random walks on the Koch network family, based 
on which we calculate analytically the average sending time (AST) defined as the average of MFPTs from a 
hub node to all other nodes, excluding the hub itself. The obtained closed-form expression displays that in 
large networks the AST grows with network order as A" In N, which is larger than the linear scaling of ART 
to the hub from other nodes. On the other hand, we also address the case with the information sender 
distributed uniformly among the Koch networks, and derive analytically the global mean first-passage 
time, namely, the average of MFPTs between all couples of nodes, the leading scaling of which is identical 
to that of AST. From the obtained results, we present that although hub nodes are more efficient for 
receiving information than other nodes, they display a qualitatively similar speed for sending information 

■ as non-hub nodes. Moreover, we show that that AST from a starting point (sender) to all possible targets is 
not sensitively affected by the sender's location. The present findings are helpful for better understanding 
random walks performed on scale- free small- world networks. 
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PACS. 05.40.Fb Random walks and Levy flights - 89. 75. He Networks and genealogical trees - 05.60.Cd 
Classical transport - 05.10.-a Computational methods in statistical physics and nonlinear dynamics 

1 Introduction cesses and underlying architecture of complex networks, 

i.e., unravel how deeply the structural features of net- 
In recent ten years, as a powerful mathematic tool, as well works affect dynamical processes occurring on them. It 
as a paradigmatic model in the intense research of com- has been shown that the power-law degree distribution of 
plex systems, complex networks have attracted a surge of scale-free networks fundamentally influence almost all dy- 
interest from the scientific community p E I M I] • Most en- namical processes taking place on them, such as disease 
deavors in the initial few years were devoted to unveil the spreading [5] , percolation [5] , games [TU1ITT] , synchroniza- 
$_i . nontrivial topological properties of real systems Q]|2]. A tion [T2], to name a few. 
^ lot of empirical studies unraveled that a large variety of 

real-life networks display simultaneously small-world ef- In addition to above-mentioned dynamical processes, 
feet [5] and scale-free behavior characterized by a power- scale-free structure also strongly affects the efficiency for 
law degree distribution [5J. These two important discover- random walks with an immobile trap fixed at a hub node 
ies have radically altered our understanding for structural with the highest degree mm[T3[THl[T71[IB] ' K W&S SUr " 
aspects of complex networked systems. prisingly found that the average receiving time (ART), 
After making substantial progress in characterizing the Le -> the average of mean first-passage time (MFPT) for a 
complexity of real systems, the focus has shifted to dy- random walker to a given target hub node, averaged over 
namical processes defined on them 0, with the aim to a11 source P oints in scale-free small- world networks, be- 
uncover the intrinsic relationship between dynamical pro- haves sublinearly or linearly with the network order (viz., 
the number of all nodes) . Here the MFPT from site utoii 
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Fig. 1. (Color online) Iterative construction method for the 
Koch networks. 




u to first reach v [T9ll20l . Since the random walker can be 
looked upon as an information messenger [2011 21] , the low 
ART to the hub node means that as information receivers 
nodes with large degree are efficient in receiving informa- 
tion. However, any node in a network can also be treated 
as an information sender. Then, interesting questions are 
raised naturally: What is the scaling of the average send- 
ing time (AST), defined as the average of MFPTs from a 
hub node to any other node, chosen uniformly in a scale- 
free network? Is it still as efficient as the case that the hub 
is regarded as a receiver? Does the location of information 
sender affect the scaling of AST? Despite the significance 
of the questions, they still remain unclear limited by the 
difficulty for determining MFPT from a hub node to some 
other nodes [22] . 

In this paper, we study analytically random walks on 
the class of Koch networks with scale-free behavior and 
small- world effect [2"3"1[M] , which is a fundamental process 
gaining considerable recent attention 25 , 26 ,27,28,29,30, 
|3"TU3"2"II3"3"] . We first investigate a particular random walk, 
starting from a hub node with highest degree to send in- 
formation to all other nodes, exclusive the hub itself. We 
derive exactly the AST from the hub to another node, av- 
eraged over all nodes in the Koch networks. The obtained 
explicit formula displays that in large networks with N 
nodes, the AST grows asymptotically with TV as N\aN, 
which in sharp contrast to the linear dependence of the 
ART from all nodes to the hub [2"3] . 

In the second part of this work, based on the connec- 
tion between random walks and electrical networks, we 
determine analytically the global mean first-passage time 
(GMFPT), defined as the average of MFPTs over all node 
pairs. We present that the GMFPT is also asymptotic to 
N In N. Since the GMFPT can be looked upon as the aver- 
age of ASTs with the sender distributed uniformly among 
all nodes, we conclude that neither the structure inhomo- 
geneity nor the position of starting points has an essential 
effect on the scaling of AST in Koch networks. Thus, the 
N\nN behavior of AST from a particular sender is a rep- 
resentative property of the Koch networks, which is in 
comparison with the trapping problem, where the scaling 
of ART to a trap (information receiver) depends on the 
location of the trap [T7] . 



Fig. 2. (Color online) A network corresponding to the case of 
m = 2. 

2 Koch networks and their structural 
properties 

The family of Koch networks controlled by a positive in- 
teger parameter m are translated from the famous Koch 
fractals [31] and can be built in an iterative way [2"51l24| . 
Denote by K mit the Koch network family after t iterations. 
Then, the Koch networks can be created in the following 
way: Initially (i = 0), K m o consists of three nodes form- 
ing a triangle. For t > 1, if m ,t is obtained from K m ^ t -\ by 
adding m groups of nodes for each of the three nodes of ev- 
ery existing triangle in K m ^-\- Each node group includes 
two nodes, both of which and their "mother" node are 
linked to each other constituting a new triangle. In other 
words, in order to get K m}t from K m j-i, one can substi- 
tute a connected cluster on the right-hand side (rhs) of 
arrow in Fig. [T]for each triangle in K m> t-\- Figure [5] illus- 
trates a Koch network for the case of m = 2 after several 
iterations. 

By construction, we can obtain with ease some quan- 
tities that will be very useful for deriving the basic quan- 
tity we are concerned in this paper. It is obvious that 
the number of triangles L&(t) present at iteration t is 
L&(t) = (3m + 1)*, and the number of nodes generated 
at iteration t is L v (t) = 6mL&(t — 1) = 6m (3m + l) t_1 . 
Then, the numbers of edges and nodes in networks K mit 
are 

E t = 3L A {t) = 3(3m + 1)' (1) 



and 



N t = V L v (ti) = 2 (3m + !)* + ! 



ti=0 



(2) 



respectively. 

Denote by fcj(i) the degree of a node i at iteration 
t that entered the networks at iteration (step) t$ (ij > 
0). Then, ki(U) — 2. Denote by L&(i,t) the number of 
triangles passing by node i at step t. According to the 
network generation algorithm, each triangle passing node 
i at a given step will lead to m new triangles involving i at 
next time step. Hence, L&{i,t) — (m + 1) L&(i,t — 1) = 
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(m + 1)* In addition, the relation fe,(t) = 2L&(i,t) 
holds. Then fcj(t) = 2(m + l) t_ * ; that indicates 



fci(t) = (m + l)fei(t-l). 



(3) 



Note that in K m j the initial three nodes created at iter- 
ation have the highest degree 2(m + 1)*. We call these 
nodes hub nodes and label by 1 one of the hub nodes, 
while label the other two hubs by 2 and 3, respectively. 

The Koch networks exhibit some classic characteris- 
tics of real-life systems [23,24 . They are scale- free with 
their degree distribution P{k) following a power-law form 



P(k) ~ k 7 , where 7 is equal to 1 



ln(3ro+l) 



belonging 



ln(m+l) 

to the interval [2, 3]. They display small- world effect with 
a small average path length (APL) and a large clustering 
coefficient. Their APL exhibits a logarithmic scaling with 
network order N t . 



3 Basic properties of Random walks on Koch 
networks 

After introducing the Koch networks K m ^ and their topo- 
logical features, we proceed to study standard random 
walks [2D] running on K mit . At each step the walker, lo- 
cated on a given node, moves uniformly to any of its near- 
est neighbors. Our main aim is to find the AST from one 
of the three hub nodes (e.g., node 1) to another node av- 
eraged over all target nodes except the hub node itself. To 
achieve this goal, we provide some essential properties for 
random walks on the Koch networks. 

3.1 Evolutionary rule for mean first-passage time 

We fist establish the scaling relation governing the evo- 
lution of MFPT between an arbitrary pair of two nodes, 
using the approach based on underlying backward equa- 
tions [22H3^1[5rJ] . Let (t) express the MFPT of the walker 
in networks if m ,t, starting from node i to visit node j for 
the first time. Because of the particular construction the 
Koch networks, the exact relation governing F^ (t+1) and 
Fij (t) can be given. 

Consider an arbitrary node i in the Koch networks 
K m ,t after t iterations. Equation @ indicates that upon 
growth of the networks from generation t to t+1, the 
degree fcj(i) of node i grows by to times, i.e., it increases 
from ki(t) to (to + l)k l {t). Denote by X the MFPT from 
node i to any of its fcj old neighbors belonging to K m j, 
and denote by Y MFPT for a walker starting from any of 
the mki new neighbors of i created at iteration t+1 to one 
of its ki old neighboring nodes previously existing before 
iteration 4+1. Then the following simultaneous equations 
hold: 



x = — W + -22,(1 + y). 

m+1 m+1 V /' 

y = ±(i + A) + ±(i + y), 



(4) 



which result in X = 3m + 
grow from generation t to t - 



1. Thus, when the networks 
- 1, the MFPT from any node 



i (i £ K m . t ) to any node j (j G 
average 3m times, namely, 



K, 



increases on 



Fij(t+1) = (3m+l)F y (t). 



(5) 



This scaling is a basic property of random walks on the 
Koch networks, which is very useful for deriving our main 
result. 



3.2 Scaling relation and expression for average return 
time 

Let Ri(t) denote the expected time for a walker in net- 
works K m j originating from node i to return to the start- 
ing point i for the first time, named mean return time 
(MRT) in the following text. By definition, we have 



h(t) 



(6) 



where J?j is the set of neighbors of node i, which belong 
to K m ,t- 

On the other hand, for K m t+i, 



Ri(t+1) 



1 



1 



TO + I 



1 h(t) 



3 en. 



(7) 

which can be elaborated as follows. The first term on the 
rhs of Eq. ([7J describes the process where the walker moves 
from node i to its new neighbors and back. Since among 
all i's neighbors belonging to K m t +i, ^rj- of them are 
new, such a process happens with a probability of 
and takes three time steps. The second term on the rhs of 
Eq. ([7]) accounts for the process in which the walker steps 
from i to one of the old neighbors j previously existing in 
K mj t and back; this process occurs with the complimen- 
tary probability = 1 - ^xx- Using Eqs. © and © 
to simplify Eq. ([7]), we can obtain the following relation 



Ri(t+1) 



3m + 1 
m + 1 



Ri(t). 



(8) 



We next determine the MRT for an arbitrary newly 
born node in K myt that is generated at iteration t. Let %' be 
a new neighbor of an old node i existing in K m>t -\, which 
is created at iteration t. Note that when i' was generated, 
another new node i" appeared at the same time and is 
linked to i and i' . Let A express the MRT of a walker 
starting off from i in networks K m j without ever visiting 
i' and i". Then we have the following relations 

Ri> (*) = ~ [1 + Fw (*)] + ~ [1 + F VH , (t)} , (9) 
F m >{t) = \ x l + i[l + F^(t)], (10) 



and 

F iV {t) 



1 ' + r^TT [1 + F V , V (t)] + k -%-^ [A + F H , (t)} . 



h{t) ki{ty 



h(t) 



(11) 
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The three terms on the rhs of Eq. (jTTj) can be under- 
stood based on the following three processes: with proba- 
bility jg-fWr , the walker gets from node i to i' in one time 

step; with probability , the walker reaches node i" in 
one time step then takes time -Fj"j'(t) to visit i'; and with 
the remaining probability fc 'ff t "7 2 , the walker selects uni- 
formly a neighbor node except i' and i" and spends on 
average time A in returning to i then takes time Fui (t) to 
arrive at node i'. 

In order to close Eqs. @ and ([IT]), we write the MRT 
of node i as: 



1 



h(t) 



x 3 + 



1 



x 3 + 



h{t) 



h(t) 



x A. 



fl2) 



The first (second) on the rhs of Eq. (|12D describes the 
process that the walker steps from i to i' (i") and back, 
which occurs with probability and needs three time 
steps. The explanation of the third term is analogous to 
that of Eq. ([TTj). 

Eliminating the three intermediate quantities Fai(t), 
Fim'{t), and A, we have 



Rl ,{t)= k M Ri[t) . 



(13) 



Combining Eqs. (0) and (fT3|) and considering fcj(i) = 2(m+ 
ly-ti j ea( j ^ following closed-form expression 



= 3(3m + l)*. 



(14) 



Note that Eq. (|14D can also be obtained from the Kac 
formula [37,38.., which states that the MRT for a node is 
in fact the inverse probability to find a particle at this 
node in the final equilibrium state of the random-walk 
process. 

Equation (TUT) does not depend the degrees of the old 
nodes, to which the new nodes i' is connected, which 
means that all the simultaneously emerging new nodes 
have identical MRT. Since all nodes born at the same time 
step have identical degree, this is obvious from the Kac 
formula: for any node with degree k in K mtt , its MRT is 
which is consistent with Eq. (| 14|) and in turn implies 
that Eq. (HH is right. 



4 Average sending time from a hub node to 
another node selected uniformly in the Koch 
networks 

In this section, we investigate the AST from a hub node 
to another node distributed uniformly in the Koch net- 
works. We focus on the case that the starting point is the 
hub node 1. Notice that, due to the symmetry, the starting 
position can be also node 2 or node 3, which does not have 
any influence on the AST. In what follows, we will show 
that the particular selection of the starting point makes it 
possible to derive analytically the relevant quantity, i.e., 
AST from node 1 to all other nodes. Let Ti(t) express the 



MFPT of node i in K m> t, which is the expected time for a 
walker starting from node 1 to first hit node i. The aver- 
age of MFPT Ti(t) over all target nodes in K m ,t is AST, 
presented by (T) t , the explicit determination of whose so- 
lution is a main goal of the following text. 

For the sake of convenient description for calculating 
(T) t , we use A t to denote the set of nodes in K m ^ t , and 
use A t to present the set of nodes created at generation t. 
Thus, we have A t — A t UA t -i. By definition, the quantity 
concerned (T) t can be defined as 



(T)t 



1 



A* — 1 



Ttot(t), 



(15) 



where T tot (t) is the sum of MFPTs for all nodes starting 
from the hub node 1, i.e., 



(16) 



ieA t 



Thus, the problem of determining (T) t is reduced to 
finding T to t(t). Since all nodes in K m j belong to either 
Z\f_i or A t , T io t(t) can be written as the sum of the two 
following terms: 

T tot (t)= T At)+ E T ^)- 

j'eA t jeAt-i 

Using the relation provided by Eq. ([SJ, Eq. pT|) can be 

recast as 

TtotW = y, T i'W + ( 3m + 1 ) Tto *(* - x ) ■ ( 18 ) 

j'eA t 

Hence, to calculate T to t(t), one only need to evaluate 
the first term on the rhs of Eq. ([TS]) , which accounts for 
the sum of the MFPTs from node 1 to all newly generated 
nodes at step t. Since before visiting node j' for a walker 
starting from node 1, it must first arrive at node j (an old 
neighbor of j') that previously existed at step t — 1, then 
Tji(t) can be written as: 



T j ,(t)=T j (t)+F jj ,(t) 



(19) 



Next we will show that Fjj/(t) can be expressed in 
terms of the quantity Rj>(t) that has been determined 
in preceding section. Note that when j' was born, it was 
linked to node j and a simultaneously emerging node j" 
that was also connected to j, then we have the following 
useful relations: 



and 



R r {t)= l -[l + F nl {t)]+ l -[l + F r , r {t)]. 



F j »f(t) = l + l[l + F jj ,(t)} 



Plugging Eq. (|2T) into Eq. O leads to 



(20) 
(21) 

(22) 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 



5 



Inserting the obtained result for Fjj> (t) given in Eq. (|22l) 
into Eq. (fl~9]) . we obtain 



Inserting Eq. (gSJ) into Eq. (fig]) , we can solve Eq. (1181) 
inductively to yield 



T f (t)=T j {t) + -R j ,(t) 



(23) 



With the result given by Eq. ([23]) . the first term on the 
rhs of Eq. ([18]), denoted by T$(t), can be expressed as 



!&>(«) = Y T r {t)= Y, [T j (t) + -R j ,(t)-2 



T tot (t) = -(3m+l)*- 1 [(12r7rf + 12™ + 2)(3m + l) i 



j'eA t 



j'eAt 



(24) 

Since for any node j created at step tj that belongs to 
At-x, there are L&(j,t — 1) = (m + 1)* _ *j _1 triangles 
passing by j, each of which will lead to 2m new nodes 
connecting j at step t, then using Eqs. (fl~4l and Eq. (f2"4"|) . 

the sum (t) can be rewritten as 

T$(t) = Y 2miA(3 > *-l)T J -(t) + 



(JV t -JV t _i) 



- x 3(3m + l)' - 2 



^ 2m(m + l)*-^- 1 T J (t) 



(JV t -N t -i) 



- x 3(3771 + 1)* - 2 



(25) 



The second term on the rhs of Eq. ([25]) is easy to 
compute. So, we only need to work out the first term 
on the rhs of Eq. (|25p . represented by T sum (i), namely, 
T 8Um (t) = Ej 6 /i«-i 2m ( m + l)* - *' -1 ^). Evidently, we 
have the following recursive relation 

^sum (t) 

= (3m+l)(m+l) Y 2m(m+l) t - t i- 2 T j (t-l) 
+2m(3m+l) Y ^H* -1 ) 

= (3m + l)(m + l)T sum (f - 1) + 2m(3m + l)T^(i - 1) . 

(26) 

On the other hand, Eq. (|25t can be rewritten as 

T$(t) = T 8Um (t) + (N t - JV t _ x ) (4(3m + 1)* - 2) . (27) 

Considering the initial conditions 7^(1) = 96m 2 + 20m 
and T sum (l) = 24m 2 + 8m, we can solve recursively the 
simultaneous equations (|26l) and (12~T1) to obtain 

T sum (i) = 8m (3m + l)*" 1 +3m(2t-l)(3m + l) 2(t " 1) , 

(28) 

and 

T^t(t) = 4m(3m + l)* _2 [6(2mf + 2m + l)(3m + if 

-3m -1]. (29) 



-3mt — 3m + ll 



(30) 



Inserting Eq. ((30]) into Eq. (]T5l) . we obtain the explicit 
expression for the AST (T) t : 



(T)t 



3(3m + 1) 



[(12mi+12m+2)(3m + l)*-3mt-3m+l] 

( 31 ) 

We continue to show how to express the key quan- 
tity (T)t in terms of the network order Nt, in order to 
obtain the relation between these two quantities. Recall- 
ing Eq. ©, we have (3m + 1)' = (N t ~ l)/2 and t = 
[ln(N t - 1) - In 2]/ ln(3m + 1). Thus, Eq. J3T]) can be fur- 
ther expressed as a function of Nt as 



(T)t = 



N t - 1 /12m[ln(A f - 1) - In 2] 
3(3m + 1) V ln(3m + 1) 
2 /3m[ln(iV t -l)-ln2] 



3(3m + 1) 



ln(3m + 1) 



f 12m + 2 
3m — 1 



(32) 



Thus, for large networks, 



(T) t 



4m 



(3m + l)ln(3m+ 1) 



(7V t -l)ln(A t -l), (33) 



showing that the AST grows with increasing order Nt as 
Nt In Nt. This leading asymptotic N t \nN t dependence of 
AST on the network order is in contrast with the linear 
scaling of receiving efficiency on network order for a re- 
ceiver located at the same hub node receiving information 
sent from all other different nodes [25] . 

It is known that the exponent 7 of degree distribu- 
tion for a scale-free network characterizes the inhomo- 
geneity of the network, which often strongly affects the 
dynamical processes running on the network [3j4,7 . As 
shown in section [2J the exponent in the Koch networks is 
7=1 + ^(m+i) , implying that parameter m controls the 
extent of heterogeneous structure of the Koch networks: 
the larger the value of m, the more heterogeneous the net- 
works. However, as shown in Eq. (|3"3"]) . although for differ- 
ent m the AST of whole family of Koch networks is quan- 
titatively different, it exhibits the same scaling behavior 
despite the distinct extent of structure inhomogeneity of 
the networks corresponding to m. 



5 Global mean first-passage time for the 
broadcaster uniformly distributed among all 
nodes 

In the previous section, we have presented that the AST 
from a most connected node to another node, averaged 
over all possible target points, exhibits a linear dependence 
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with network order by a logarithmic correction. However, 
for this case, the information sender is placed on a largest 
node. Then a question arises naturally whether this scaling 
is representative. Another interesting issue is whether the 
diffusion speed still follows the same behavior when the 
sender is located on other nodes. In the following text, 
we will study the case that the information sender is uni- 
formly distributed among all nodes in the networks, in 
order to explore how deeply the position of the sender 
affect the scaling of transportation efficiency. 



5.1 Exact solution to global mean first-passage time 

In this case, we are concerned in a new quantity called 
global mean first-passage time (GMFPT), which is the 
average of mean first-passage times over all pairs of nodes 
in the networks. Concretely, the GMFPT in K m j, repre- 
sented by (F)t, is defined as 



(F)t 



in which 



1 



Ftotjt) 

N t (N t - 1) N t (N t - 1) 



N t N t 



n EE ( 34 ) 



N t N t 



(35) 



3 = 1 ?=i. 

in 



is the sum of MFPTs between all pairs of nodes. Note that 
the definition of GMFPT involves a double average: The 
first one is over all the walkers to a given target (receiver) 
j, the second one is over a uniform distribution of target 
nodes among all nodes in K m j. 

It should be noticed that the above method used for 
computing (T) t is not applicable to (F) t , so we must re- 
sort to an alternative approach. Fortunately, the peculiar 
construction of the Koch networks and the link [39,40 
between effective resistance and the MFPTs for random 
walks allow to calculate analytically GMFPT (T) t . We 
view K mi t as resistor networks [41] by considering each 
edge to be a unit resistor. Let Rij(t) be the effective resis- 
tance between two nodes i and j in the electrical networks 
obtained from K m t - Then, according to the relation be- 
tween MFPTs and effective resistance [391140] . we have 

F ij (t)+F ji (t) = 2E t R ij (t). (36) 

Therefore, Eq. ([55} can be rewritten as 

N t N t 

-Aot 

(*) = ^EE i2 *i(*)- ( 37 ) 

j= i «=i 

Thus, if one knows how to determine the effective resis- 
tance, then we have a method to find (F) t . Then, the 
question of determining (F) t is reduced to computing the 
total resistance i?tot(i) between all pairs of nodes in the 
resistor networks: 



Aft N t 



«tot(*)=EE J2 *i(*)- 



(38) 



According to the structure of the Koch networks, it 
is obvious that the effective resistance between any two 
nodes is exactly § times the usual shortest-path distance 
between the corresponding nodes, i.e., 



(39) 



where dij (t) is the shortest distance between nodes i and 
j in K m ,t- Equation (I39[) can be interpreted as follows. 
By construction, the Koch networks consist of triangles; 
moreover, no edge lies in more than one triangle. Then, 
for any couple of nodes i and j in K m j, the shortest path 
between them is unique. It is easy to see that the effective 
resistance between two nodes directly connected by an 
edge in the shortest path of i and j is |, which is in fact 
equal to the effective resistance between two nodes of a 
triangle. And the Rij(t) can be regarded as the sum of 
effective resistance of dij(t) conductors in series, each of 
which has a effective resistance of |. 

Then, to obtain (F) t , we need only to calculate the 
total of shortest distances between all node pairs, denoted 
by Aot(*), namely 



Nt 



It is then obvious to have 



Rtot(t) = -Aat(t). 



(40) 



(41) 



Hence, all that is left to find (F) t is to evaluate D tot (t). 

According to our previous result [24;, we can easily 
obtain the closed-form expression for D tot (t): 



Aot(*) 
2(3m + I)*" 1 



"3m + 5 + (24mi + 24m + 4) (3m + 1)'] . 

(42) 



Combining above-obtained results, we arrive at the ex- 
plicit solution to (F) t : 



(F)t 
2 I 



3N t (N t -l) 
2(3m + I)*" 1 



6(3m + I)' + 3 



Et Aot(t) 

[3m + 5 + (24mt + 24m + 4) (3m + 1)'] , 



(43) 



which can be expressed in terms of network order Nt as 
1 N t -1 



(F)t = 



3(3m + 1) N t 
12mln(7V t - 1) - 12m In 2 
ln(3m + 1) 



3m + 5 + (N t - 1) 
12m + 2 



(44) 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 



7 



Equation (|44D uncovers the exact dependence relation 
of GMFPT on network order N t and parameter m. For 
large systems, i.e., Nt — > oo, we have following expression 
for the leading term of (T) t : 



(F)t 



Am 



(3m+ l)ln(3m + 1 



-(JVt-l)KNt-l), (45) 



which is in consistent with the general result given in |17j . 
Thus, similar to the behavior of AST obtained in the pre- 
vious section, in the large limit of t, the GMFPT grows 
with network order Nt as NlnNt, which is independent 
of m and thus shows that the structure heterogeneity of 
the networks has no substantial impact on the scaling of 
GMFPT. The sameness for the leading behavior between 
{T)t and (F) t implies that the N t In N t scaling of (T) t from 
a hub node to all other nodes is a representative feature 
for information sending in the Koch networks. 

The N t In N t behavior found for both the AST and 
GMFPT can be understood from the following heuris- 
tic explanations. The couples of nodes farthest apart (be- 
tween each other and from the hub due to its centrality) 
provide the leading contribute for the related MFPTs [42 . 
On the other hand, for the ART related to the trapping 
problem with the trap fixed on a hub node, since the hub 
is relatively easy to reach for most nodes, the ART is rela- 
tively small and contributes little to GMFPT, see also [TB] . 

Note that if the information sender is positioned at 
an arbitrary non-hub node in networks K m j. The AST 
from the sender to all other nodes also follows the scal- 
ing Nt In Nt ■ Because in most of this case, the information 
must be first delivered to a hub node in a time at most 
proportional to network order Nt [23], then the piece of 
information proceeds to be sent, until it reaches the re- 
ceiver after an average transmit time Nt In Nt as shown in 
the previous section. To confirm this, we have computed 
analytically the AST for the sender located at new neigh- 
bor of hub node 1 created at step t, and obtained the same 
expression as Eqs. (J33J and (|4"S"j) . 



6 Conclusions 

We have studied random walks on the Koch network fam- 
ily, exhibiting synchronously scale-free and small-world 
behaviors. We first concentrated on a specific case for ran- 
dom walks from a hub node to all other nodes, and ob- 
tained explicitly the formula for AST from this most con- 
nected node to different target nodes, which varies with 
network order N as TV In TV, larger than the ART from 
all other nodes to the hub. Then we continued to derive 
the GMFPT between two arbitrary nodes averaged over 
all node couples in the Koch networks, which can be re- 
garded as the average of MFPTs from a uniformly-selected 
starting point to all other nodes in the networks. We pre- 
sented that in the limit of large network order N, the 
GMFPT also scales approximatively with TV as TV In TV. 
This identity of scalings between the AST and GMFPT 
indicates that the ability (efficiency) of hub nodes sending 
information is the same as that the average efficiency of 



all nodes in the Koch networks, showing that the sending 
efficiency measured by AST is not sensitively influenced 
by the position of information sender and the structural 
heterogeneity of the networks. Finally, it should be men- 
tioned that we only studied a particular family of scale- 
free networks, whether the conclusion also holds for other 
scale-free networks, even general networks, needs further 
investigation in the future. 
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